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Abstract 

Smoking and COPD are associated with decreased mucociliary clearance, and healtliy smol<ers have shorter cilia in the large 
airway than nonsmokers. We hypothesized that changes in cilia length are consistent throughout the airway, and we further 
hypothesized that smokers with COPD have shorter cilia than healthy smokers. Because intraflagellar transport (IFT) is the 
process by which cilia of normal length are produced and maintained, and alterations in IFT lead to short cilia in model 
organisms, we also hypothesized that smoking induces changes in the expression of IFT-related genes in the airway 
epithelium of smokers and smokers with COPD. To assess these hypotheses, airway epithelium was obtained via 
bronchoscopic brushing. Cilia length was assessed by measuring 100 cilia (10 cilia on each of 10 cells) per subject and 
Affymetrix microarrays were used to evaluate IFT gene expression in nonsmokers and healthy smokers in 2 independent 
data sets from large and small airway as well as in COPD smokers in a data set from the small airway. In the large and small 
airway epithelium, cilia were significantly shorter in healthy smokers than nonsmokers, and significantly shorter in COPD 
smokers than in both healthy smokers and nonsmokers. The gene expression data confirmed that a set of 8 IFT genes were 
down-regulated in smokers in both data sets; however, no differences were seen in COPD smokers compared to healthy 
smokers. These results support the concept that loss of cilia length contributes to defective mucociliary clearance in COPD, 
and that smoking-induced changes in expression of IFT genes may be one mechanism of abnormally short cilia in smokers. 
Strategies to normalize cilia length may be an important avenue for novel COPD therapies. 
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Introduction 

The lung is constantly exposed to inhaled particles, pathogens, 
irritants and toxins. In the normal lung, airway mucus provides a 
barrier to the airway epithelium and traps inhaled particles, which 
are removed from the airway via mucociliary clearance, the 
process by which cilia beat in a coordinated fashion to sweep 
mucus up out of the airways [1-3]. Impaired mucociliary 
clearance leads to a cycle of excessive airway mucus, recurrent 
pulmonary infection and worsening obstructive airway disease 
[4,5]. In chronic obstructive pulmonary disease (COPD), muco- 
ciliary clearance is dysfunctional for a variety of reasons, which 
may include alterations in mucus composition and generation, and 
adverse effects of cigarette smoke on cilia structure and function 
[5,6]. Prior work in our laboratory has shown that in the large 
airway epitheKum (LAE, 2"''-5''^ order airways), cigarette smoking 
in healthy individuals is associated with shorter cilia, with an 
observed difference in length that is likely to have an effect on 
mucociliary clearance [7]. 



How cilia length is controlled is not well understood. A number 
of theoretical mechanistic explanations have been proposed. These 
include the molecular ruler model, in which a protein with a 
physical length matching the cilia length controls cilia length [8,9]; 
the limited precursor model, in which cilia length is limited by 
limited quantities of necessary precursor molecules [9,10]; the 
cumulated strain model, in which cilia length is controlled by 
binding energy changes induced by changes in conformation with 
increasing length [9,11]; a feedback control model, involving 
length sensing and signal transduction [9,12]; and the balance 
point model, which is based on the knowledge that ciliary 
disassembly at the tip is ongoing at steady-state and appears to be 
length-independent, and which states that there is specific length 
at which this disassembly process is balanced with the rate of 
assembly, dependent on the rate of intraflagellar transport of 
components to the cilia tip [9,13]. 

In addition to global mechanistic explanations, a number of 
specific genes have been implicated in length control. Some of 
these genes have been discovered through genetic screening of 
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Chlamydomonas flagellar length mutants, including LF2, a cyclin- 
dependent kinase family member, and LF4, a mitogen-activated 
protein kinase family member [14]. In our previous report of 

shorter cilia in the large airway epithelium in liealthy smokers, we 
found decreased expression of DNAH5, — 10, and — 11, as well as 
DYNC2H1, ezrin, and ODF2 in healthy smokers as compared to 
nonsmokers [7]. Consistent with the balance point model of length 
control, other investigators have specifically assayed proteins 
involved in intraflagellar transport, including those related to both 
the anterograde kinesin motor and the retrograde dynein 2 motor, 
as well as the intraflagellar transport genes IFT27, IFr46, IFr52, 
IFT70, IFT88 and IFT172 [13,15-23], showing in model 
organisms that these genes are necessary for the production 
and/ or maintenance of cilia of normal length. 

We hypothesized that the changes in cilia length observed in the 
large airway of smokers arc global throughout the airway, and that 
COPD smokers have shorter cilia than smokers without evidence 
of COPD. In exploring how smoking might impact ciUa length, 
based on published observations of the central role of IFT in cilia 
length, we hypothesized that smoking causes decreased expression 
of intraflagellar transport genes in the human airway epithelium 
and that this reduced expression is associated with the shorter cilia 
observed in smokers. We addressed these hypotheses by assessing 
cilia length in the large and small airway epithelium of 
nonsmokers, healthy smokers, and smokers with COPD, and by 
quantifying expression of 40 IFT genes in the airway epithelium of 
these groups. 

Methods 

Ethics Statement 

All individuals were evaluated and samples collected in the 
Weill Cornell NIH Clinical and Translational Science Center and 
Department of Genetic Medicine Clinical Research Facflity under 
cUnical protocols approved by the Weill Cornell Medical CoUege 
and New York/Presbyterian Hospital Institutional Review Boards 
(IRB) according to local and national IRB guidelines. AU subjects 
gave their informed written consent prior to any clinical 
evaluations or procedures. 

Study Population 

Healthy nonsmokers, healthy smokers, and smokers with 
COPD were recruited from the general New York City population 
by placing advertisements in local newspapers and electronic 
bulletin boards. Individuals were determined to be phenotypicaUy 
normal or to have COPD based on history, physical exam, 
complete blood count, coagulation studies, liver function tests, 
urine studies, chest X-ray, EKG, and pulmonary function tests. 
COPD was defined and staged according to the GOLD criteria 
[24]. Smoking status was confirmed with urine nicotine and 
cotinine levels. Full inclusion/exclusion criteria are detailed in 
Text SI. For COPD subjects, exacerbation frequency. Modified 
Medical Research Council dyspnea score [25], and St. George's 
Respiratory Questionnaire scores [26] were assessed. High 
resolution chest CT was performed on a subset of subjects for 
quantitative analysis of emphysema (see Text SI). This study is 
registered under the ClinicalTrials.gov identifiers NCT00224185 
and NCT00224198. 

Sample Collection and Preparation 

AU subjects underwent bronchoscopy for the acquisition of 
airway epithelial ceUs using standard methods [27-29]. FuU details 
are available in Text SI. 



Cilia Length Analysis 

CUia length was assessed in the LAE and SAE of healthy 
nonsmokers, healthy smokers and COPD smokers in a random- 
ized, blinded fashion. One slide was evaluated per subject. Each 
slide was visualized via bright field microscopy under 60x 
magnification using a Nikon Microphot-SE (Nikon, Melville, 
NY) and images of 10 ciliated cells per subject were taken using an 
Olympus DP-70 color CCD camera (Olympus, Center VaUey, 
PA). Using Image J software (NIH, 2011), images were contrast 
enhanced and 1 0 cilia on each of 1 0 cells were measured for each 
subject. A mean cilia length was calculated for each individual. 
The average of aU the means for each subject within a phenotype 
was taken for the group mean. Standard error of the mean and 
coefficient of variation were calculated for each individual and 
each phenotypic group. Further details regarding the method for 
determining mean cilia length and statistical analysis are provided 
in Text SI, Figure SI and Table SI. 

IVlicroarray Data Analysis 

Microarray data were processed using the MASS algorithm 
(Aflymetrix Microarray Suite Version 5 software) and GeneSpring 
software. FuU details are available in Text SI. The raw microarray 
data are pubUcaUy avaUable at the Gene Expression Omnibus 
(GEO) site (http://www.ncbi.nlm.nih.gov/geo/) under the acces- 
sion number GSE43939. Changes in the expression of significant 
genes were confirmed with TaqMan PGR; fuU detaUs are available 
in Text SI. 

Results 

Study Population 

Ciha length was assessed in 228 airway epithelium samples, 
including 120 LAE samples (n = 25 healthy nonsmokers, n = 25 
healthy smokers, and n = 70 COPD smokers [n = 34 GOLD I, 
n = 29 GOLD II, n = 7 GOLD III]) and 108 SAE samples (n = 20 
healthy nonsmokers, n = 32 healthy smokers, and n = 56 COPD 
smokers [n= 18 GOLD I, n = 29 GOLD II, n = 8GOLD III, 
n = 1 GOLD IV]) (Table 1). Gene expression was assessed in the 
LAE of nonsmokers (n = 21) and healthy smokers (n = 31) and the 
SAE of an independent group of nonsmokers (n = 28) and healthy 
smokers (n = 69) as weU as in the SAE of smokers with COPD 
(n = 59; n = 38GOLD I, n = 20GOLD II, n=l GOLD III) 
(Table 1). Both the healthy smokers and the smokers with COPD 
were generaUy older than the nonsmokers (gene expression set 
LAE p>0.1, SAE p<0.01; cilia length set LAE p<0.001, SAE 
p<0.003). Among subjects with LAE samples, healthy smokers 
and COPD smokers reported similar pack-yr of smoking (p>0.05); 
however, among subjects with SAE samples, smokers with COPD 
had higher reported pack-yr than healthy smokers (p<0.02 in the 
cUia length set and p<0.001 in the gene expression set). In the 
gene expression data set, among the LAE samples, there were no 
significant differences in the cell differentials between healthy 
smokers and nonsmokers (aU p>0.1). Among the SAE samples, 
smokers had significantiy fewer cfliated and basal ceUs and more 
undifferentiated ceUs compared to nonsmokers (aU p<0.03); 
smokers with COPD had fewer ciliated ceUs than healthy smokers 
(p<0.02). In the cilia length data set, among LAE samples, COPD 
smokers had significantly more epithelial cells and fewer inflam- 
matory cells than nonsmokers (both p<0.02), though all pheno- 
types had >98% epithehal ceUs. Healthy smokers and COPD 
smokers had fewer cUiated ceUs than nonsmokers, and COPD 
smokers had more secretory cells than nonsmokers and healthy 
smokers and more undifferentiated ceUs than nonsmokers (aU 
p<0.05). Among the SAE samples, COPD smokers had fewer 
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ciliated cells than nonsmokers and healthy smokers, more 
undifferentiated cells than nonsmokers, and more basal cells than 
healthy smokers (all p<0.04). In both the LAE and SAE groups 
smokers with COPD were slightly taller; this difference reached 
statistical significance (p<0.03) when compared to nonsmokers in 
the LAE group only. COPD smokers with GOLD II and III 
disease had significandy more emphysema on HRCT than those 
with GOLD I disease (5.1% vs 2.8%, p<0.02). 

Effect of Smoking and COPD on Cilia Length in the Large 
and Small Airway Epithelium 

Cilia length was measured as described in Methods for 
nonsmokers, healthy smokers and COPD smokers (Figure S2). 
The method of measuring 100 cUia on one slide for each 
individual was used after determining that mean cUia length was 
not significandy different when 100 or 500 cUia were measured, 
and that intra-slide variability for one subject was less significant 
than intra-individual variability (see Text SI and Table SI). 
Greater variability in cilia length for a given individual was 
observed among COPD smokers and healthy smokers than among 
nonsmokers (see Text 81 and Figure S3). The mean length for 
each phenotype was calculated as a mean of the individual means. 
In both the LAE and SAE, healthy smokers had shorter cUia than 
nonsmokers, and cilia were shorter in COPD smokers than in 
nonsmokers or healthy smokers (Figure 1). In the LAE, mean cUia 
length was 7% shorter in healthy smokers as compared to 
nonsmokers (7.09 vs 7.63 |J,m), 12% shorter in COPD smokers as 
compared to healthy smokers (6.16 vs 7.09 |a,m), and 19% shorter 
in COPD smokers as compared to nonsmokers. In the SAE, mean 
cilia length was 9% shorter in healthy smokers as compared to 
nonsmokers (6.49 vs 7.15 [un), 6% shorter in COPD smokers as 
compared to healthy smokers (6.05 vs 6.49 |a.m), and 15% shorter 
in COPD smokers as compared to nonsmokers. Mean cilia length 
was shorter in the SAE than in the LAE with significant differences 
observed in both nonsmokers (7.15 vs 7.63 (xm) and healthy 
smokers (6.49 vs 7.09 Jim), with a similar trend in COPD smokers 
(6.05 vs 6A6 \lm). 

To exclude an effect of age or sex on the findings, subsets of 
nonsmokers, healthy smokers and COPD smokers were chosen 
that were matched for age and sex. In the LAE, each group 
consisted of 4 women and 9 men and mean age was 47 ±7 for 
nonsmokers, 48 ±6 for healthy smokers, and 48 ±5 for COPD 
smokers (p>0.8). Mean cilia length was 7.8 [wa for nonsmokers, 
7.4 |im for healthy smokers, and 6.2 ^im for COPD smokers 
(p<0.001 for COPD vs nonsmoker, p<0.002 for COPD vs 
smoker, but p>0.4 for smoker vs nonsmoker). In the SAE, each 
group consisted of 6 women and 9 men and mean age was 
45±9 for nonsmokers, 47±6 for healthy smokers, and 47±6 for 
COPD smokers (p>0.5). Mean cflia length was 7.1 |xm for 
nonsmokers, 6.6 |J.m for healthy smokers, and 5.9 |J,m for COPD 
smokers (p<0.001 for COPD vs nonsmoker, p<0.03 for COPD vs 
smoker, but p>0.1 for smoker vs nonsmoker). 

Because the COPD smokers in the SAE group of the study 
reported significandy more pack-yr of smoking than did the 
healthy smokers, and the COPD smokers in the LAE group had a 
nonsignificant trend toward more pack-yr of smoking than healthy 
smokc-rs, a subset of COPD smokers in each group was randomly 
selected that was matched for pack-yr with the healthy smokers. In 
the LAE, the matched subset of COPD smokers (n = 25) had a 
mean age of 49 ±5 compared to 49 ±5 for healthy smokers (n = 25, 
p>0.9) and mean pack-yr 33±11 vs 32±10 for healthy smokers 
(p>0.9). Mean cUia length for this subset of COPD smokers was 
6.19 [im vs 7.09 |im for healthy smokers (p<0.001). Similarly, in 
the SAE, the matched subset of COPD smokers (n = 32) had a 



mean age of 50 ±7 yr, compared to 48±5 yr for the healthy 
smokers (n = 32, p>0.1) and a mean smoking history of 
29± 10 pack-yr, compared to a mean smoking history of 
29± 1 1 pack-yr for the healthy smokers (p>0.9). In this subset 
analysis, the mean cUia length remained significandy shorter in the 
COPD smokers: 6.01 (im vs 6.49 [im in healthy smokers (p<0.02). 

Clinical and Demographic Correlations with Cilia Length 

Correlations between mean cilia length and demographic and 
lung function parameters were assessed using univariate linear 
regression analysis (for continuous variables) and ANOVA (for 
categorical variables). In the LAE analysis, in the total population, 
there were significant correlations between cilia length and age, 
ethnicity, FEVl, FVC, FEVl/FVC, and DLCO (all p<0.01). 
When these factors were incorporated into a multivariate analysis, 
only FEVl, FVC, and FEVl/FVC remained as significant factors 
(all p<0.02). This analysis was repeated excluding nonsmokers 
and analyzing smokers (healthy and COPD) only. In the 
univariate analysis, significant factors were age, sex, ethnicity, 
FEVl, FEVl/FVC, and DLCO; in the multivariate analysis, only 
the FEVl remained as a factor significandy correlated with mean 
cilia length (p<0.0001). In univariate analysis of FEVl vs mean 
cUia length including all subjects, the r^ value was 0.25 (Figure 2A). 
COPD smokers with GOLD II-III disease had slightly shorter 
mean cilia length than those with GOLD I disease, but this 
difference was not significant (6.1 |im vs 6.2 [im, p>0.4). 

Analysis was repeated in the same manner for the SAE. In the 
total population, significant correlations existed between cilia 
length and age, ethnicity, height, FEVl, FVC, FEVl/FVC, and 
DLCO (all p<0.03). In multivariate analysis, height and FEVl/ 
FVC remained significant (both p<0.005). When only smokers 
were considered, significant correlations were found for age, sex, 
ethnicity, height, FEVl, F\'C, and FEVl/FVC, but none of these 
correlations remained significant in multivariate analysis. In 
univariate analysis of FEVl/FVC in all subjects, the r^ value 
was 0.19 (Figure 2B). Consistent with the relationships observed 
between cilia length and spirometric variables, subjects with 
GOLD II-IV COPD had significantly shorter cilia than those with 
GOLD I disease (5.9 [im vs 6.4 |j,m, p<0.01). 

Effect of Smoking on Intraflagellar Transport Genes 

To evaluate the hypothesis that smoking-induced changes in 
cilia length may be due, in part, to changes in the expression of 
genes encoding proteins relevant to intraflagellar transport, a list of 
40 known, human IFT-related genes was compiled and expression 
of these genes evaluated in the LAE and SAE [30]. All 40 known 
IFT-related genes were expressed in both the LAE and SAE data 
sets (Table S2). Analysis of expression in smokers vs nonsmokers in 
the LAE revealed 1 3 probe sets representing 1 1 genes with 
expression significandy modified in smokers. Of these 1 1 genes, 
10 were down-regulated and 1 (IFT80) was up-regulated in 
smokers. The down-regulated genes included a component of the 
anterograde IFT motor, and 4 genes involved in the IFT complex 
B (anterograde transport), including TRAF3IP1 gFT54), IFT57, 
IFT172, and CLUAPl. DYNC2H1, a component of the 
retrograde transport motor, was down-regulated in smokers, as 
were 2 components of the IFT complex A (retrograde transport), 
IFT43 and WDR19 (IFT 144). Two components of dre BBSome, 
an IFT complex accessory, were also down-regulated, BBSS and 
BBS9 (Figure 3A). 

Expression analysis in the SAE revealed similar results. Sixteen 
probe sets representing 14 genes were significantly modified in 
smokers; of diese, 3 were up-regulated (HSPBl 1/IFT25, TTC26, 
and BBIPI/BBIPIO) and 11 down-regulated. Among the antero- 
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grade transport category, down-regulated genes were KIF3A and 
KIF17, as well as TRAF3IP1 (IFT54),' IFT57, IFT172, and 
CLUAPl. Down-regulated genes also included the retrograde 
transport gene DYNC2H1, as well as another retrograde motor 
component, WDR34, and IFT43, WDR35 (IFT121), and 
WDR19 (IFT144; Figure 3B). 

Combining these two analyses revealed a total of 8 IFT-related 
genes for which expression was similarly down-regulated by 
smoking in the independent data sets of LAE and SAE. These 
included 5 anterograde transport genes, KIF3A, TRAF3IP1/ 
IFT54, IFT57, IFT172, and CLUAPl; and 3 retrograde transport 
genes, DYNC2H1, IFT43, and WDR19/IFT144. TaqMan PCR 
was carried out for significant genes on samples from the LAE and 
confirmed the changes observed by microarray (Table S3) with a 
significant p value for 5 of these genes. 
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Effect of COPD Status on Intraflagellar Transport Genes 

To assess whether the difference in mean cilia length in COPD 
smokers compared to healthy smokers could be due to further 
down-regulation of IFT genes in this group, we assessed IFT gene 
expression in COPD smokers compared to smokers without 
COPD. No significant differences in gene expression were found 
in COPD smokers vs healthy smokers for the 40 IFT genes 
assessed. We also assessed for differences in expression of these 
genes between GOLD I and GOLD II-III subjects and found no 
differences. 

Correlation of Gene Expression with Cilia Length 

Finally, we assessed whether the expression of the 8 IFT genes 
with smoking-induced down-regulation as confirmed in two 
independent data sets correlated with cilia length. This was done 
using the subset of subjects with SAE samples for whom both cilia 
length and gene expression were available (n = 83; there were no 
LAE subjects who had both gene expression and cilia length data). 
Genes were assessed first using univariate analysis. The results 
indicated that for 5 of the 8 genes, there was a modest but 
significant correlation between expression and cilia length: KIF3A, 
TRAF3IP1, IFT57, IFT43, and WDR19 (all p<0.05; Figure 4A- 
E). These 5 genes were then incorporated into a multivariate 
analysis and this model also significandy correlated with cilia 
length (p<0.005, r2 = 0.19). 

Discussion 

Disorders of cilia are typically thought of as inherited disorders, 
the most prominent cHiopathy being primary ciliary dyskinesia, a 
clinical syndrome composed of pulmonary disease, male infertility, 
and, in 50% of patients, situs inversus [31,32]. COPD can be 
thought of as an acquired, rather than inherited, cUiopathy in 
which cigarette smoke induces dysfunction in mucociliary clear- 
ance by a variety of mechanisms related both to abnormal mucus 
and mucin expression [5,6], as well as to abnormal cilia structure 
and function [6,33,34]. This mucociliary dysfunction can be 
devastating, as mucus stasis leads to a vicious cycle of recurrent 
respiratory infections and worsening airflow obstruction [4,5], 
Proposed mechanisms of dysfunctional cilia in smokers include 
decreased numbers of ciliated cells, absent cilia structural 
components, changes in cilia beat frequency and coordination of 
ciliary beating, and altered mucus hydration status [34—37]. We 
have reported that healthy smokers have shorter cilia in the large 
airway than do nonsmokers [7] . In the current study we show that 
smokers with COPD have even shorter cilia than healthy smokers 
and that these changes are present in the distal as well as in the 
proximal airway, and we present data suggesting that smoking 



PLOS ONE I www.plosone.org 



5 



January 2014 | Volume 9 | Issue 1 | e85453 



IFT Genes and Cilia Length 




Figure 1. Cilia lengtKi in tlie large and small airway of nonsmolters, healthy smokers, and COPD smokers. IVIean cilia length in 
micrometers is displayed on the ordinate. The abscissa displays large airway on the left and small airway on the right, with healthy nonsmokers 
represented in white, healthy smokers in dark gray, and smokers with COPD in black. Error bars show the standard error of the mean for each 
phenotype. Within the large airway, n = 1 20 (healthy nonsmokers n = 25, healthy smokers n = 25, smokers with COPD n = 70). Within the small airway, 
n = 108 (healthy nonsmokers n = 20, healthy smokers n = 32, smokers with COPD n = 56). There are significant differences in mean cilia length 
between all phenotypes in both the large and small airway (p<0.05). There are significant differences in mean cilia length between airway locations 
in both healthy nonsmokers and healthy smokers (p<0.05). 
doi:1 0.1 371 /journal.pone.0085453.g001 



represses the expression of genes involved in intraflagellar 
transport, which may be at least partially causative of the short 
cUia, and that expression of several IFT-related genes correlates 
with cilia length. 

Very few studies have addressed the effect of smoking on airway 
epithelial cilia length. In 1957, Chang [38] evaluated respiratory 
epithelium obtained at autopsy from the lungs of smokers with 
unknown lung function (n=71) and nonsmokers (n = 34) and 
observed that smokers had shorter airway cilia than nonsmokers. 
Serafmi and Michaelson [39] observed that cilia length decreases 
with increasing airway order in subjects without lung disease (but 
unknown smoking status) and described shorter cilia in a single 
subject with COPD. Tamashiro et al [40] found that exposure to 
cigarette smoke extract inhibited cUiogenesis and resulted in 
shorter cilia in cultured mouse nasal epithelium. Most recently, we 
reported that cilia in the large airway epithelium of clinically 
healthy smokers were significandy shorter than those of healthy 
nonsmokers [7]. 

To understand potential mechanisms for the shorter cUia 
observed in smokers and COPD smokers, we focus(-d on 
intraflagellar transport (IFT), the process by which molecular 
motors transport and assemble tubulin at the distal ends of cilia 



[30]. IFT is required both for the growth and maintenance of 
normal cilia, as cilia maintenance is a dynamic process with 
ongoing turnover of tubulin at the distal tip necessitating transport 
of cUia proteins to the distal tip (anterograde IFT) and back to the 
cell body (retrograde IFT) [30] . Anterograde transport is driven by 
the kinesin-2 motor and utilizes proteins of IFT complex B. 
Retrograde transport is dri\ (;n by the multiprotein cytoplasmic 
dynein 2 motor and utilizes proteins of IFT complex. 

A number of proteins related to IFT have been shown in model 
systems to directly affect cilia length. Disruption of the anterograde 
kinesin motor results in short or absent cilia in a number of 
organisms [13,30] as do alterations in the retrograde cytoplasmic 
dynein motor [15,30]. Homologs of IFT27, IFT46, IFT52, IFT70, 
IFT88 and IFT172, all relevant to anterograde transport, as well 
as the BBSome component BBS9, have been demonstrated in 
model organisms to be necessary for generating and/ or maintain- 
ing cUia of normal length [16-23,41]. Less is known about genetic 
control of cUia length in humans. Interestingly, however, members 
of a consanguineous family with short-rib Polydactyly syndrome 
were found to have mutations in the dynein heavy chain 
DYNC2H1 resulting in decreased mRNA expression, and 
cultured chondrocytes from these individuals demonstrated short 
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Figure 2. Correlation between mean cilia length and lung function. Each triangle represents one subject, with nonsmokers represented by 
green triangles, healthy smokers by yellow triangles, and COPD smokers by red triangles. A. Correlation of FEV1 with mean cilia length in large airway 
epithelium. The ordinate represents FEV1 (% predicted) and the abscissa represents mean cilia length in the LAE. B. Correlation of FEV1/FVC with 
mean cilia length in the small airway epithelium. The ordinate represents FEV1/FVC (% observed) and the abscissa represents mean cilia length in the 
SAE. 
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cilia and cells with absent cilia [42] . In our previous study of cilia 
length in the large airway of smokers, DYNC2H1 was down- 
regulated in the healthy smokers [7] . Further, members of a family 
with Sensenbrenner syndrome, a ciliopathy characterized by 
abnormal bony development as well as frequent extraskeletal 
abnormalities, were found to have a mutation in WDR19 
(encoding IFT 144); cultured fibroblasts from a Sensenbrenner 
patient also displayed reduced cilia length compared to cilia from 
fibroblasts from control cell lines [43] . 

To expand on our observation of shorter cUia in the LAE of 
healthy smokers, we sought to address two questions: (1) how does 
the cilia length of smokers with COPD compare to that of 
nonsmokers and healthy smokers; and (2) are these changes also 
seen in the small airway, the earliest site of abnormality in the 
development of COPD? The data indicate that the cUia of COPD 
smokers are 6 to 12% shorter than the cUia of smokers without 
COPD and 15 to 19% shorter than the cilia of healthy 
nonsmokers. This reduction in length would be expected to 
significantly adversely affect mucociliary clearance, both by 
reducing the ability of the cUia tips to make contact with the 
mucus layer as well as by reducing the force generated on the 
forward stroke of the cilia, and provides an additional mechanism 
by which mucociliary clearance is impaired in smokers with 
COPD. Our results also show that smokers and COPD smokers 
have shorter cUia in both the large (3"^'^-4* order) and small ( 1 0*- 
12*) order airways. Statistically significant correlations were 
observed between mean cilia length and physiologic markers of 
airflow obstruction, FEVl in the LAE analysis and FEVl/FVC in 
the SAE analysis, suggesting that cUia length may play a role in the 
pathogenesis of disease in smokers. In the SAE analysis, a 
significant correlation was also found for height. The significance 
of this is uncertain, and the degree of correlation was small 
(r^ = 0.05); this is unlikely to have influenced our results, as COPD 
subjects in our study were on average slightly taller than 
nonsmokers. 

To assess potential mechanisms for the observed differences in 
cUia length, we hypothesized that smoking down-regulates 
expression of IFT genes in the airway epithelium and used 



microarray technology to assay the expression of 40 human IFT 
genes representing both anterograde IFT (the kinesin-2 motor and 
IFT complex B) and retrograde IFT (the cytoplasmic dynein 2 
motor, IFT complex A) as well as the BBSome in the human 
airway epithelium in smokers and nonsmokers in two independent 
data sets representing large and small airway epithelium. The 
results indicated 8 IFT genes with expression modulated by 
smoking as confirmed in the two independent data sets. Among 
genes involved in anterograde transport, 5 were significantiy 
down-regulated in association with smoking in both data sets. 
These included both genes for which the literature supports a role 
in maintenance of cilia length (KIF3A [13] and IFT172 [23]), and 
genes for which no specific role in cUia length maintenance has 
been elucidated (IFT54, IFT57, and CLUAPl). Similarly, among 
the 3 down-regulated retrograde transport genes, there is data 
supporting the role of DYNC2H1 and IFT 144 in cilia length 
[15,30,42,43], and we have previously reported on decreased 
DYNC2H1 expression in healthy smokers with short cilia [7], but 
no specific role in cilia length maintenance has been reported for 
IFT43. In the context that smoking, in general, both up- and 
down-regulates large numbers of genes in the airway epithelium 
[27-29,44—46], with some categories of genes, such as oxidant- 
related genes, being strongly up-regulated by smoking, the finding 
of down-regulation of IFT genes supports the concept that there is 
a specific effect of smoking on IFT. The mechanism by which 
smoking affects both anterograde and retrograde IFT genes is an 
important area for future study. One potential mechanism is 
smoking-induced suppression of a common cUia-related transcrip- 
tion factor, such as FOXJ 1 . However, other mechanisms, such as 
epigenetic modification, may be involved, and the present study 
does not answer this question. 

Interestingly, despite the finding of shorter cUia in COPD 
smokers as compared to healthy smokers, even when matched for 
pack-yr, no differences in IFT expression were observed between 
these groups. The present study does not explain why this is the 
case. Given the complexity of cilia assembly and maintenance, and 
the many effects of cigarette smoke on the airway epithelium, this 
finding may suggest that there are additional mechanisms aside 



PLOS ONE I www.plosone.org 



7 



January 2014 | Volume 9 | Issue 1 | e85453 



IFT Genes and Cilia Length 



o 
o 



o 
ji. 



o 



o 

bo o 



Relative gene expression 

K ^ rv 



O 

b 



O 



O 



70 

> 



o 

c 
> 

■0 



O 

-< 



o - 













o 








CJ 


Tl 




H 




-(^ 




OJ 


§ 




D 




73 


H 


w 


H 




O 




N) 




CD 






O 










CD 


CD- 



CD . 

w 



03 
00 



□ A 
p 
O 

■o 

□ A 
° 
s 

J p 

p 




— » o 

■D 

□ s 

o 
■a 

□ V 
p 
u> 
o 



3 

0) 



a) 
»< 

<D 
3- 

E" 
3 



(D 

I 

■o 



—J o 



-1 

9 
8 

□ A 
g 



— ' o 



W 

2 3 
a 



E" 
3 



Figure 3. Expression of intraflagellar transport (IFT)-related genes significantly modified by smoldng in tKie airway epithelium. Gene 
expression was assessed in nonsmolcers and healthy smokers by microarray analysis in two independent data sets. The ordinate represents relative 
gene expression and specific genes are displayed on the abscissa. Nonsmokers are represented by light bars and healthy smokers by dark gray bars. 
Error bars represent standard error of the mean and p values are calculated using the Benjamini-Hochberg correction. Significant p values are 
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identical, the probe set with highest fold-change is displayed. A. Expression in the large airway epithelium (nonsmokers n = 21 and healthy smokers 
n = 31). B. Expression in the small airway epithelium (nonsmokers n = 28 and healthy smokers n = 69). 
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from alterations in IFT contributing to cilia length control in 
smokers. Alternatively, given the variability in cUia length observed 
among individuals, it may suggest that those smokers with the 
shortest cilia are at increased risk for the development of COPD. 

The role of cUia in cells is diverse and is not limited to a 
clearance function. Emerging data suggests a role for motile cilia 
in cell signaling and IFT proteins may play a direct role in several 
signaling pathways [47]. An intriguing area for future study is the 
effect of smoking on the signaling function of airway cUia. 

These results support the concept that loss of normal cilia length 
may contribute importantly to the defective mucocUiarv' ck-arancc 
seen in COPD. The fact that abnormally short cUia are also found 
in the small airways of healthy smokers and COPD smokers is 
consistent with the obstruction and dysfunctional mucociliary 



clearance seen in the small airways in COPD [4,5]. In conjunction 
with the findings of short ciUa in large and small airways of healthy 
and COPD smokers, we have found that 8 genes involved in IFT 
are down-regulated in the airway epithelium in association with 
smoking, and the expression of 5 of these genes correlates with cilia 
length. These results do not prove causation; such studies, for 
instance evaluating the effect of knocking down expression of IFT 
genes in human cells, would be very difficult given the multiple 
genes involved in IFT with probable overlapping function. 
However, while these results must be interpreted with caution, 
the inclusion in this list of 4 genes with known or suggested roles in 
maintenance of normal cUia length adds strength to the findings. It 
is likely that the effect of smoking on IFT is complex and involves 
the suppression of a number of key IFT genes, which may be one 
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mechanism of abnormally short cilia in healthy and diseased 
smokers. Further study of the process by which smoking results in 
shorter cilia in smokers and of how these shorter cilia impact 
mucociliary clearance should lead to novel therapeutic avenues 
targeted at the dysfunctional mucociliary clearance that is a key 
feature of COPD. 
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